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We analyze the infrared behavior of effective iV-point interactions between order 
parameter fluctuations for nematic and other quantum critical electron systems with 
a scalar order parameter in two dimensions. The interactions exhibit a singular 



1} . momentum and energy dependence and thus cannot be represented by local vertices. 
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\ They diverge for all N > 4 in a collinear low-energy limit, where energy variables 

scale to zero faster than momenta, and momenta become increasingly collinear. The 
i ■ degree of divergence is not reduced by any cancellations and renders all iV-point 

interactions marginal. A truncation of the order parameter action at quartic or any 

o. 

other finite order is therefore not justified. The same conclusion can be drawn for the 
£N| ■ effective action describing fermions coupled to a U(l) gauge field in two dimensions. 
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I. INTRODUCTION 

Numerous interacting electron systems undergo a quantum phase transition^ between 
ground states with different symmetries, which can be tuned by a non-thermal control 
parameter such as doping, pressure, or a magnetic field. In the vicinity of a continuous tran- 
sition electronic excitations are strongly scattered by critical order parameter fluctuations, 
such that Fermi liquid theory breaks down.-^ Quantum critical fluctuations near a quantum 
critical point (QCP) are therefore frequently invoked as a mechanism for non-Fermi liquid 
behavior in strongly correlated electron compounds. 

Quantum criticality in metallic electron systems is traditionally described by an effective 
order parameter theory which was pioneered by Hertz^ and extended to finite tempera- 
tures by Millis.- In that approach an order parameter field cf) is introduced via a Hubbard- 
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Stratonovich decoupling of the electron-electron interaction, and the electronic variables are 
subsequently integrated out. The resulting effective action S[4>] for the order parameter is 
truncated at quartic order and analyzed by standard scaling techniques. 

However, several studies revealed that the Hertz-Millis approach may fail, especially in 
low-dimensional systems.- 1 ^ Since electronic excitations in a metal are gapless, integrating 
out the electrons may lead to singular interactions between the order parameter fluctuations 
which cannot be approximated by a local quartic term. The nature of the problem and es- 
sential aspects of its solution were presented first for disordered ferromagnets by Kirkpatrick 
and Belitz. 7 For clean ferromagnets, Belitz et al.- showed that Hertz-Millis theory breaks 
down, and no continuous quantum phase transition can exist, in any dimension d < 3; the 
transition is generically of first order. 9 The Hertz-Millis approach was also shown to be in- 
valid for the quantum antiferromagnetic transition in two dimensions.— ~— In that case a 
continuous transition survives, but the QCP becomes non-Gaussian. 

In this article we analyze the validity of the Hertz-Millis approach to quantum critical- 
ity for two-dimensional systems with singular forward scattering of electrons, in particular 
systems exhibiting a quantum phase transition driven by forward scattering in the charge 
channel. The most prominent such transition is the electronic nematic, in which an ori- 
entation symmetry is spontaneously broken, while translation and spin-rotation invariance 
remain unaffected.— The problem of quantum critical points with singular forward scatter- 
ing is formally similar to the problem of non-relativistic fermions coupled to a U(l) gauge 
field, which was studied intensively in the 1990s.— ^ 

Perturbation theory for the electronic self-energy at the nematic QCP yields a non- Fermi 
liquid contribution proportional to tu 2 ^ 3 already at the lowest order in a loop expansion.— 1 ^ 
The same behavior was found earlier for fermions coupled to a U(l) gauge field.— It was 
commonly believed that the power-law with an exponent | is not modified by higher order 
contributions. Furthermore, calculations in the gauge field context suggested that the sim- 
ple form of the (bosonic) fluctuation propagator obtained in lowest order (RPA) remains 
unaffected by higher order terms. 19 In a remarkable recent paper Metlitski and Sachdev 20 
formulated a scaling theory of the nematic QCP and related problems, treating the electrons 
and order parameter fluctuations on equal footing. In a renormalization group calculation 
they found a logarithmic divergence at three-loop order pointing at a correction of the cu 2 ^ 3 
law for the electronic self-energy. However, no qualitative correction was found for the flue- 
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tuation propagator, up to three-loop order.— & This is in stark contrast to the case of an 
antiferromagnetic QCP in two dimensions, where the fluctuation propagator is substantially 
renormalized compared to the RPA form.— ^ A clarification of the properties of the nematic 
QCP beyond three-loop order is still lacking. 

The robustness of the fluctuation propagator at the nematic QCP seems to indicate that 
interactions of the order parameter fluctuations are irrelevant such that the QCP is Gaussian, 
in agreement with the expectations from Hertz-Millis theory. It is therefore worthwhile to 
analyze the interaction terms in the effective action S[<f>] obtained after integrating out the 
electrons. The A-point interactions are given by fermionic loops with N vertices. To obtain 
the scaling behavior of such loops is non-trivial, because the most naive power-counting is 
easily invalidated by cancellations.— 1 ^ In this paper we compute the exact scaling behavior 
of the A-point interactions for the nematic QCP and related systems. We find that the 
interactions are marginal and non-local for all A > 3. Hence, replacing them by a local (/> 4 
interaction is not justified. 

The paper is structured as follows. In Sec. II we introduce the effective actions to be 
analyzed, and we define the A-point loops describing the interaction terms. In Sec. Ill 
we explain the special role of fluctuations with collinear momenta, which motivates the 
definition of the collinear low-energy scaling limit. Sees. IV- VI are dedicated to the analysis 
of the A-point loops. After reviewing exact formulae from the literature (Sec. IV), we derive 
explicit expressions for iV-point loops in the collinear low-energy scaling limit. Using the 
scaling behavior of these loops, we perform the power counting of A-point order parameter 
interactions in Sec. VII. We finally conclude in Sec. VIII. 



We consider an interacting Fermi system which undergoes a continuous quantum phase 
transition with a scalar order parameter of the form 



where c£(k) and c CT (k) are the usual fermionic creation and annihilation operators. For 
a charge nematio^ the form factor d CT (k) is spin symmetric and has a k-dependence with 
G?-wave symmetry, such as d a (k) = cosk x — cosk y . A spin-antisymmetric form factor may 



II. EFFECTIVE ACTION AND N-POINT LOOPS 




(1) 
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describe an Ising ferromagnet or an Ising spin nematic. 

Decoupling the fermionic interaction by introducing an order parameter field <p y i a a 
Hubbard-Stratonovich transformation, and integrating out the fermionic variables,- one ob- 
tains an effective action 

S[<j>] = \f g-'^qm-q) 
00 (-1)^ f 

N = 2 Jqi,...,q N 

where g > is the fermionic coupling constant, and 




Here and in the following we use 3-vectors collecting imaginary frequency and two- 
dimensional momentum variables, for example k = (/c ,k), and j k as a short-hand notation 
for J ^(f^y?- The variables pj and qj are related by 

1j = P3+1 ~Pj for 3 = !> • • • > N ~ 1 

qN = pi-pn- (4) 

Note that qi + ■ • • + q^ = due to energy and momentum conservation. The bare propagator 
has the form Go(k) = [iho — e(k) + , where e(k) is the dispersion relation of the non- 
interacting particles. Hd,N(qi, ■ ■ ■ , qN) can be represented graphically as a fermion loop with 
N lines corresponding to Go and N vertices with form factors d a (k), as shown in Fig. 1. For 
spin-antisymmetric form factors, II^at vanishes for odd N. 

For fermions coupled to a (7(1) gauge field, integrating out the fermions leads to a similar 
effective action S[<f>], where <fi is the transverse component of the gauge field (in Coulomb 
gauge). The bosonic TV-point functions are then essentially given by a loop with transverse 
current vertices 

N 

n t ,jv(gi, ...,q N ) = 2 / ] J [q j± ■ v(k - Pj - qj/2) G (k - Pj )} , (5) 
Jk j=i 

where v(k) = Ve(k), and qj_ is the unit vector obtained by rotating q = q/|q| by ir/2, that 
is, qj_ = (—q y ,q x )- In addition there are contributions from the "diamagnetic" term (of the 
form (f) 2 ipil}) in the underlying fermionic action, which are however less singular, since they 
involve a smaller number of propagators (less than TV). 
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FIG. 1: Graphical representation of Il^jV with momentum variables as in Eq. ([3]). 

The iV-point contribution to S[<f>] in Eq. (j2]) is symmetric under any permutation of 
qi, . . . , qN- Hence, one can replace Hd,N by the symmetrized iV-point loop 

n d,jv(9i> • • • , Qn) = -^j- ^2 U d,N(qpi, ■ qpN) , (6) 

p 

where the sum collects all permutations P of 1,...,N. Non-symmetric contributions to 
n^jv do not contibute to the integral in Eq. (j2]). For the gauge field problem one defines 
nf N analogously. Substantial cancellations may occur in the sum over permutations.— 1 ^ 



III. COLLINEAR LOW-ENERGY LIMIT 

The quadratic part of S[4>] is determined by a constant and the 2-point loop or "bubble" 
U d (q) = U d , 2 (q, ~q) = J2 [ G (k - q/2)G (k + q/2) . (7) 

For small q and small |?o|/|c[|) it h as an expansion of the form^i 

Mq) = -N d + X ,q 2 + 7d M + . . . , ( 8 ) 

where N d = J j0n d 2 a (k)5[e(k) — /x] is a weighted density of states, and \d and 7d are two 
other constants. At the quantum critical point one has g^ 1 — N d — such that the quadratic 
part of the action vanishes for q — > and <Zo/|q| ~~ * 0- F° r the gauge problem, the constant 
N t from the static limit of n t (g) cancels generically against a tadpole contribution,—^ such 
that the theory is always critical. In both cases, the Gaussian part of the action has thus 
the asymptotic form 

SoM = \f (W + 7^y) <Kq)<K-q) , (9) 
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corresponding to a bare propagator 

A»(g)= * M - (10) 

Do(q) diverges in the limit q — > and go/|q| — > 0. The two terms in the denominator of 
Do(q) vanish at the same pace for q — y if go oc |q| 3 . Therefore, the bare dynamical scaling 
exponent is z = 3. To assess the size of the interaction terms in S[<f>] one thus has to study 
the iV-point loops in a low-energy limit with qjo oc |q/| 3 . Naively one would expect that this 
corresponds to the static limit, where qjo — y before q^ — y 0. Hertz and KlenirJS showed 
that an iV-point loop converges to the (N — 2)-th derivative of the density of states with 
respect to the Fermi energy in the static limit. In our case, with a form factor d a (k), their 
result generalizes to 

i im ^o n ^ (91 ' • • • ' «») = jN^iy. aj?^ 22 J j^f $ ( k ) ~ A ■ ( n ) 

Except for special cases where the chemical potential lies at a van Hove singularity, this 
expression is finite. Note that the right hand side of Eq. ffTTj) is independent of qi, . . . , q^ and 
hence already symmetrized. Approximating the bosonic iV-point interactions by finite local 
interactions thus seems adequate. Standard power counting then implies that all interactions 
with iV > 4 are increasingly (with higher N) irrelevant. Hence, the Hertz- Millis truncation 
seems justified for QCPs with singular forward scattering, even in two dimensions. The 
static limit of the 3-point loop and all other iV-point loops with odd iV as given by Eq. ffTTj) 
vanishes even in the case of a charge nematic, due to the antisymmetry of d a (k) under 7r/2 
rotations of k. 

One arrives at a similar conclusion for the gauge field problem. In that case the static limit 
of Ht,N is also generically finite. Odd iV-point interactions vanish due to the antisymmetry 
of v(k) and even iV-point interactions appear to be irrelevant for any iV > 4. A Gaussian 
fixed point thus seems natural. 

However, there is a flaw in the above argument. Eq. fTTTl) has been derived by setting 
qjo = before the momenta q 7 - tend to zero. It is not guaranteed that this captures the 
low-energy limit qj — y and <2jo/|qj| - > in general. Indeed, a simple estimate indicates 
that the iV-point loop is of order g/o/lq/l^" 1 for small non-collinear momenta cy and small 
9jo/|Qj|-~ Although this behavior is increasingly singular for larger N, the corresponding 
order parameter interactions remain irrelevant, since the singularity is not strong enough.— 
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However, an even stronger singularity is obtained in a special low-energy limit in which 
the momenta qi,...,qjv become collinear. The crucial role of coupled fluctuations with 
collinear momenta was highlighted very clearly by Metlitski and Sachdev.— In perturbative 
one-loop calculations of the fermionic self-energy £(kp, u) at a certain point kp on the Fermi 
surface, it was found already some time ago that the dominant contributions involve only 
fermionic states in the momentum region near kp and — kp , with momentum transfers q 
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for an 



almost tangential to the Fermi surface in those points. See, for example, Ref. 
early calculation in the gauge field context, and Ref. |l7j for a corresponding calculation at 
the nematic QCP. This remains true for higher order contributions,— so that all fermionic 
momenta are close to kp and — kp and momentum transfers are almost tangential to the 
Fermi surface in these points, which implies that they are mutually almost collinear. 

Choosing a coordinate system in momentum space in such a way that the normal vector 
to the Fermi surface at kp points in x-direction, the proper scaling limit describing the 
low-energy behavior is given by ko i-> X 3 ko, k x h-> X 2 k x , and k y h-> Xk y with A — > 0, where 
(k x , k y ) is measured relative to kp.— >2Z>2& For the momentum and energy transfers qj this 
implies 

q j0 H> A 3 g iQ , q jx H> X 2 q jx , q jy H> Xq jy (12) 

with A —7- 0. In this collinear low-energy limit the momentum transfers q, become increas- 
ingly collinear (pointing in y-direction). The behavior of the iV-point interactions given 
by n rf A r(g 1 , . . . , q N ) and Yl t ^(qi, ■ ■ ■ , <?a0 in the collinear low-energy limit has not yet been 
studied systematically. In particular, it has not yet been analyzed whether cancellations 
suppress their value below the naive power counting estimate. To clarify this issue is the 
main purpose of our article. 

As mentioned above, the dominant contributions are due to momenta k close to those 
points ±kp on the Fermi surface at which the momentum transfers q., are tangent. In the 
definining equation (j3J) for H^n we can therefore replace the form factors d a (k) by d a (±k F ). 
Assuming d a (— kp) = d CT (kp), which is satisfied in all cases of interest, we then obtain 

U d>N (q 1 , • • • , q N ) -> d ° ( k ^) n AK<?i, ■■■An) (13) 

17 

in the collinear low-energy limit, where 

f n 

IL N (qi,...,q N ) = I N (pi,...,p N ) = J^Goik-pj) (14) 

Jk j=i 
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is the iV-point loop for spinless fermions with unit vertices. Furthermore, the dispersion 
e(k) enters only via the Fermi velocity vp and the Fermi surface curvature in ±kp. Both 
are assumed to be finite, which is the generic case. We choose units such that vf and 
the curvature radius pf are both one, and we realize these parameters by using a simple 
parabolic dispersion relation e(k) = k 2 /2 and setting kp = 1. Relating II^jv to Il y with a 
quadratic dispersion enables us to exploit exact results for 11^ which are already available 
(see below). For the gauge field problem the vertices are antisymmetric under reflections, 
since v(— k) = — v(k). Hence, II fi jv can be reduced to II jy in the collinear low-energy limit 
only for even N: 

N 

II ti jy(gi, . . . , q N ) -»■ 2 Y[ qj± ■ v(k F ) IIjv(gi, • • • , gjv) • (15) 
i=i 

For odd N, contributions from k near kp and — hp contribute with opposite sign and can 
therefore not be written in terms of Hat- However, the results obtained for even N suffice 
to show that the effective action involves non-local marginal interaction of arbitrarily high 
order. 

The symmetrized TV-point loop 

Tl s N (qi, . . . , q N ) = U N{qpi, • • • , qpn) (16) 

p 

describes the dynamical iV-point density correlations of a Fermi gas. In the following sections 
we will derive its scaling behavior in the collinear low-energy limit for arbitrary N. 

IV. EXACT FORMULAE FOR N-POINT DENSITY LOOP 

Our analysis of the scaling behavior of n^v and is based on exact expressions derived 
by Feldman et al.— and their elaboration by Neumayr and Metzner.— They are valid for 
a parabolic dispersion relation. In this section we summarize these expressions, assuming 
specifically e(k) = k 2 /2 and kp = 1. Obviously one may restore an arbitrary mass and kp 
at will. We use the parametrization In(pi, ■ ■ ■ ,Pn) with momenta pj linearly related to the 
qj, as described in Sec. II, see also Fig. 1. The following expressions are applicable only for 
non-collinear momenta. Nevertheless, they can be used to study the scaling behavior in a 
limit where they become increasingly collinear upon reducing A. 
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The iV-point loop can be expressed as a linear combination of 3-point loops with rational 
coefficients:— 



i N (p l7 . . . ,p N ) = 

l<i<j<k<N 



N 



n /-( d 



ijk\ 



h(j>i,Pj,Pk) 



where 



(17) 



ijk\ 



+ 



«(P< ~ Pi) + KPiO ~ Pvo) 



(Pfc - Pi) + KPko ~ Pm) 



det(pj - p,,p„ - Pi 
det(p j - p l; p fc - p./ 



The determinants of two momenta are defined as det(p, p') = det y^p 
For the 3-point loop, Feldman et al.— have obtained the expression 

3 



h(pi,P2,Ps) 



1 



2vrz det(p 2 -pi,p 3 -Pi] 



S * J 



where s 
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S23 = s 3 i = 1, s 2 i = s 32 = S13 = -1, and 



7ij 



18) 



(19) 



(20) 



The contour-integrals are performed along the curves 7^ = {iOy(s)|0 < s < 1}, where Wij(s) 
is the unique (generally complex) root of the quadratic equation 



(Pj-Pi) z + 2det(d-pj,p j -p i ) z + (d-pi) = s z 
satisfying the condition 

j m ~(Pjx ~ Pix)Z + d y — 

(Pjy Piy)Z d x Pix 

The (complex) two-dimensional vector d = (d x , d y ) is given by 



(21) 



(22) 



det(p 2 - Pi,P3 - Pi 



(P3-Pi) + ^(P30-Pio) 



(P2 - Pl)± + P2 P3 ■ (23) 



The integration path UJy(s) can be written explicitly as22 
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Here = %ij(s) + iyij{s) is a function of s with real and imaginary parts given by 



x 



\J\aijis)} 2 + {Pjo-Pio) 2 + 



aij {s) 



1/2 



1 



^/[^•(s)] 2 + {p j0 -p l0 ) 2 - a i:j (s) 



1 1/2 



with 



an s 



!| ,2 

tIpj-PiI + 



(Pio - Pio) 



(25) 
(26) 

(27) 



The constant is given by = x,- t j + z?/y = x^ + iy^*., where completes the index set 
{i,j} to {i,j, k} = {1,2,3}, and 



Vijk 



I Pi - Pi 



vPj - Pk) ■ [Pk - Pi 



2det(pj - P i, Pfc - Pi 
(Pfc - Pi)(PjO - Pko) ~ (Pj ~ Pk)(PkO ~ Pm) Pj ~ Pi 



(28) 
(29) 



det(pj - Pi, p k - Pi) \Pj — Pi| 

One can easily show that yij(s) increases strictly monotonically as a function of s, and 
Xij(s) increases (decreases) strictly monotonically if sgn(p j0 ~ Pio) > (sgn(pj — Pio) < 0). 
The integration path uiij(s) thus has a simple shape. We finally note the following obvious 
symmetries under exchange of i and j: 



x ji{ s ) = -Xij(s) yji(s) = yij(s) 



Xji Xij yji x/ij 



(30) 
(31) 



V. 3-POINT LOOP IN COLLINEAR LOW-ENERGY LIMIT 

We now derive the asymptotic behavior of the 3-point loop l2,{P\iP2iPz) in the collinear 
low-energy limit. To this end we substitute pj h-> X 3 Pjo, Pj X >-> X 2 Pjx, Pjy ^ ^Pj y and 
expand I3 as given by Eq. ( fl~9l) in powers of A. 

We first expand the integration path ?%(s), Eq. (124T) . For the constants Xijk and y^ one 
obtains 

= o Hpt j^~ P n ^ ~ ~ *») + ° (A) ' (32) 

zclet^pj — p^ p k — pi) 

, (Pfcy ~ gig) (PjO ~ Pfco) - (Pjy ~ Pky) (PkO ~ Pio) Pjy ~ Piy , ,,,, . 2^ / Q o^ 

det ( pj - pi , p k - pi ) \p jy - piy I 
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For the functions Xy(s) and Vij(s) one finds 



x. 



■i j 



» = sgn(pj - Pio)s + C(A 2 



-A 



sbjo — P»ol 



2.s 



+ 0(A 5 ) 



for s > 0, and 



xy(0) = A 2 PJ0 _ ^ + Q(A 3 ) , 

3/v(0) = -^-PiJ + ^A 2 ). 

Inserting the expansion of the above auxiliary quantities into Eq. 
and imaginary parts, one obtains 

1 (Pjy ~ Pky)(Pky ~ Piy) 



(34) 
(35) 

(36) 
(37) 

and splitting the real 



Wij(0) 



and 



A 2det(p,- - PA - Pi 



+ 0(1) 



— i 



1 + (Pky ~ Piy)(PjO - Pko) ~ (Pjy ~ Pky)(PkO ~ Pio) Pjy ~ Piy + 

2 det(pj - Pi, Pfc - Pi) - Pij,| 2 



^(o)+T sgnfao : p ; o) +g(i)+^ 



(9(A) 



(38) 



(39) 



A \Pjy~Piyl 
for s > 0. 

The value of tij, Eq. ( l20l) . is given by the difference of natural logarithms at the end 
and the beginning of the integration path, plus contributions ±2iri for each crossing of the 
branch cut on the negative real axis in the complex plane. For small A one can write ty in 
a form where no case-dependent multiples of 27ti appear, namely 31 



1+i 



Uij{S) 



Uj = \n[uij(s)\ + In 
where 

The sum in Eq. (1191) can be written as 



(40) 



(41) 



E (42) 

(ij)=(l,2),(2,3),(3,l) 

Forming the difference — the first terms from Eq. (T4TJ1) cancel due to the antisymmetry 
of Uij(s) in i and j, such that 

i 



tij tji In 



l+i 



ViA s ) 

Uij{s) 



In 



1 - i 



Uij(s) 



(43) 
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This is a suitable starting point for an expansion in powers of A, since Vij(s) j 'wy(s) is of 
order A for s = 0, and of order A 3 for s — 1. Expanding the logarithm yields 



E 

71=0 



2n + 1 



! 




2n+l 


[,.y«(o)' 


2n+l ^ 













(44) 



Inserting from Eq. fT4"Tj) and expanding in powers of A, one obtains 



+ 2i 



1/^(1) 3^(0)^(0) _ 4 //,,;((') ^ (0) 



Xij{\) Xij 



X; 



x: 



3x< 



+ 0(\ A ) 



(45) 



where the first term is of order A, the second of order A 2 , and the third one of order A 3 . 

In the sum over pairs from Eq. ( 1421) . contributions of order A and A 2 to the single 

differences — tji cancel, as do many terms of order A 3 . To see this one has to insert 
expansions of the auxiliary quantities appearing in Eq. (145 j) in powers of A, sometimes beyond 
the order presented in Eqs. (152"]) - (|37|) . After a lengthy but straightforward calculation one 
obtains 3 ^ 

VijiX) , /0/x4 > 



I' 3 = 2z 



1 - Xij(l) Xij 



+ (9(A 4 ) . 



(46) 



(M)=(l,2),(2,3),(3,l) 

Expanding x^, and 1/^(1) yields 

1 ) A 3 D ijk (p j0 -p i0 ) 



where 



and 



(1) - x^ 2 + |sgn(p i0 - Pio) sgn(p jy - Pi y )F ijk 
D ijk = det(pj - p fc , Pi - pfe) , 

-^ijfc — (Pky ~~ Pjy)(Pjy ~~ Piy)(,Piy ~ Pky) ■ 



+ o{\ 4 



(47) 



(48) 



(49) 



Note that and are both invariant under cyclic permutations of and k. Inserting 
Eq. (T4T1) into Eq. (l4"6l) . and dividing by 27ri det(p2 — pi, P3 — pi), we obtain our final result 
for the collinear low-energy limit of the 3-point loop 



113(91,92,93) = h(Pl,P2,P3) 

= - E 



PjO ~ PiO 



271 (i^J^V+cyc. Dl23 + i Sgn ^'° - SgI1 ^ - ^ 123 



0(A). (50) 



13 



The 3-point loop is thus generically finite for A — > 0. The limit is real and depends on 
the ratios (p j0 -Pio)/D 123 and (p j0 - p i0 )/F 123 . Note that |F 123 | = q ly q 2 yq 3 y, while |D 123 | is 
twice the area of the triangle with corners pi, p 2 , and P3, or, equivalently, of the triangle 
obtained by attaching the vectors qi, q 2 and q 3 to each other. IT 3 (gi, g 2 , <?3) vanishes if 
frequency variables qjo are set to zero before scaling q, to zero, in agreement with the result 
of Hertz and Klenin.— It also vanishes if either sgn(jOj — p i0 ) = sgn(pj y — p iy ) for all 
or sgn(j9jo ~ Pio) — ~ s S n (Pjy ~ Piy) f° r & H (hj)- m these cases the contributions to the sum 
over cyclic permutations of (1,2,3) in Eq. (1501) cancel. 

The expression on the right hand side of Eq. ( 150|) is invariant under permutations of q%, 
q 2 , and g 3 , so that it also describes the collinear low-energy limit of the symmetrized 3-point 
loop U§(qx,q 2 ,q 3 )- 



VI. N-POINT LOOP IN COLLINEAR LOW-ENERGY LIMIT 



The reduction formula (JTTJ) relates the iV-point loop to a linear combination of 3-point 
loops. The coefficients are determined by the quantities fi V {^ k ) defined in Eq. (fTBI) . In the 
collinear low-energy limit, the latter become frequency independent and real, and they scale 
as 



+ 



2 

A 2 



i.Piy Pvy] 



(Pky ~ P 



D 



2 

'iv> D 



+ j^k 



13k 



+ 0(A 2 



(51) 



Inserting this and Eq. (|50|) for J 3 into Eq. (TT7|) . one obtains an explicit formula for the Appoint 
loop in the collinear low-energy limit. It is remarkable that Hn{qi, ■ ■ ■ , ?iv) = In(.Pi, ■ ■ ■ ,Pn) 
is a rational function of all momenta and frequencies in this limit. For A^ > 3, it diverges as 



U N oc A 2 < 3 -^ 



(52) 



for A ->• 0. 

The degree of divergence of Iljv is not reduced upon symmetrization, so that the sym- 
metrized Appoint loop also diverges as 



(53) 
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We have confirmed the absence of significant cancellations by computing the scaling behav- 
ior of n^r for various choices of iV and q±, . . . ,q^. This result is remarkable since strong 
and systematic cancellations have been shown to occur upon symmetrization when the 
limit qj — > is taken more conventionally. In particular, uniformly scaled iV-point loops 
njv(Agi, . . . , Agjv) diverge as X 2 ~ N , while their symmetrized counterparts n^(Agi, . . . , Ag^) 
remain finite.— 1 ^ The arguments establishing the cancellation of divergences in the uniform 
small-g limit do not apply in the collinear low-energy limit. 

The divergence in Eq. (|53|) is the "worst case scenario" compatible with simple power 
counting: In the collinear low-energy limit the integration measure in the definition of the 
iV-point loop scales as A 6 , while each of the N propagators diverges as A -2 , such that the 
loop may diverge as \ 6 ~ 2N (but not stronger). What we have shown is that this divergence 
is neither reduced by oscillations of the integrand under the /c-integral in Eq. f[T4"|) . nor by 
cancellations in the sum over permutations contributing to the symmetrized loops. Note that 
the power-counting does not change if the bare propagator Go(k) in the Appoint loops is 
replaced by a propagator G(k) with a self-energy proportional to fc 2 / 3 , since this interacting 
propagator also diverges as A~ 2 . 

VII. POWER COUNTING OF ORDER PARAMETER INTERACTION 

Now that we have determined the scaling behavior of the Appoint loop in the collinear 
low-energy limit, we can assess the relevance of the Appoint order parameter interactions 
in the effective action S[(/>] by using standard power counting. To see how the interaction 
terms evolve compared to the quadratic part of the action, we rescale the field <fi in such 
a way that the bare critical action So[0], Eq. (Q, remains invariant. Since the integration 
measure scales as A 6 , and the inverse bare propagator DQ 1 (q) as A 2 , we have to rescale 
the field as <p (-> A~ 4 0. The Appoint interaction terms in S[<f>] are composed of ./V energy- 
momentum integrals, a delta-function for energy-momentum conservation, the Appoint loop, 
and a product of iV fields . . . (^(Qn)- The iV-point interaction therefore scales as 




(54) 



All Appoint interactions contributing to S[<p] are thus marginal in the collinear low-energy 
scaling limit. Hence, the effective order parameter action cannot be truncated at any fi- 
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nite order. At least such a truncation is not justified by power counting. Furthermore, 
the interaction terms have singular momentum and energy dependences, which cannot be 
represented by local interactions. 

Whether the Gaussian fixed point remains stable or not depends therefore entirely on 
the behavior of fluctuation contributions. To get an idea one might compute low order 
fluctuation corrections to S/ )3 [</>] and £1,4 [</>]. Metlitski and Sachdev 2 ^ have shown that the 
bare Gaussian propagator D (q) does not receive qualitative modifications up to three- loop 
order in the coupled fermion-boson theory underlying the effective action S[4>]M^ Obtaining 
a general conclusion on the fate of the Gaussian fixed point seems difficult, however, since 
the theory has no obvious expansion parameter. 

Let us also discuss the power counting for a generalized Gaussian part of the effective 
action, of the form Sq[4>] — | J q (x ( 1 1+<e +7^)0(q , )0( — q), where e G [0, 1]. This generalization 
was introduced by Nayak and Wilczek^ 2 - for the sake of a controlled expansion in e. It was 
recently used by Mross et al.— for the purpose of defining a managable large-iV/ limit of 
the theory, where Nf is the number of fermion flavors. The case e = is related to the 
theory of electrons in a half-filled Landau level,— while e = 1 describes the nematic QCP 
and related systems as discussed above. For arbitrary e, the scaling limit f fl2|) has to be 
generalized to qjo h-> \ 2+t qjo, qj x h-> X 2 qj X , and qj y \qj y , corresponding to a dynamical 
exponent z = 2 + e. The limit is still collinear and "static" (in the sense that frequencies 
scale to zero faster than the modulus of momenta). In light of the results for e = 1 it is thus 
likely that the degree of divergence of the iV-point loops is not reduced by cancellations in 
this scaling limit for e < 1, too, such that it can be estimated by naive power counting. One 
then obtains IT^ oc \ 5+€ ~ 2N . Powercounting for the effective action then yields S7,tv[0] A 
as in the special case e = 1. The iV-point interactions are thus marginal in the collinear 
low-energy limit, for any N and e. The size of e does not matter here. 

We finally compare to the situation in three dimensions. Here finite <fr N interactions are 
even more irrelevant than in two dimensions, but one may again wonder about the special 
role of momentum transfers tangential to the Fermi surface in a certain Fermi point (and its 
antipode), which are mutually coplanar. Considering the coplanar low-energy scaling limit 
q Q 1 — y \ 3 qo, Qx l— >• A 2 g x , and q y>z h-> Xq y>z , one expects a divergence LT^ oc \ 7 ~ 2N for N > 4. 
These divergences are however not strong enough to make SV,v[0] marginal or relevant. 
Power counting yields Sj t3 [(f)\ oc A 1 / 2 , and SV,v[0] oc \ N I 2 for N > 4. The Gaussian fixed 
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point is thus clearly stable in three dimensions. 

VIII. CONCLUSION 

In summary, we have analyzed the scaling behavior of iV-point interactions between 
order parameter fluctuations at a nematic QCP and for other quantum critical electron 
systems with a scalar order parameter in two dimensions. The iV-point interactions are 
given by symmetrized fermionic loops with iV vertices. We have shown that these loops 
exhibit a singular momentum and energy dependence for all N > 3, so that they cannot be 
represented by a local interaction. For iV > 4, they diverge in the collinear low-energy limit, 
where energy variables scale to zero faster than momenta, and momenta become increasingly 
collinear. We have derived explicit expressions for the momentum and energy dependences 
in that limit. The degree of divergence is not reduced by any cancellations. From standard 
power counting one then obtains that all iV-point interactions are marginal. 

The effective action is thus dominated by interactions between fluctuations with collinear 
momenta, as noted already previously— It cannot be truncated at any finite order and none 
of the iV-point terms can be represented by a local interaction. In particular, approximating 
the interactions by a local quartic term, as in Hertz' theory, is inadequate. The same 
conclusion can be drawn for the effective action describing fermions coupled to a U(l) gauge 
field. 

Marginality of all iV-point order parameter interactions has already been obtained for 
a spin density wave QCP in two dimensions.— In that case the fluctuation propagator is 
strongly renormalized by these interactions. In particular, anomalous scaling dimensions 
appear.—^ For the nematic QCP and related theories it is presently unclear whether the 
Gaussian fixed point remains stable. Perturbative calculations have not yet revealed any 
singular renormalization of the bare Gaussian propagator.—^ 
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